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Radiative neutralino decay X2 — * Xil is studied in a Split Supersymmetric scenario, and com- 
pared with mSUGRA and MSSM. This 1-loop process has a transition amplitude which is often 
quite small, but has the advantage of providing a very clear and distinct signature: electromagnetic 
radiation plus missing energy. In Split Supersymmetry this radiative decay is in direct competition 
with the tree-level three-body decay X2 — ^ Xi //' ™d we obtain large values for the branching ratio 
B{x2 — > Xil) which can be close to unity in the region M2 ~ Mi. Furthermore, the value for the 
radiative neutralino decay branching ratio has a strong dependence on the split supersymmetric scale 
TO, which is otherwise very difficult to infer from experimental observables. 



I. INTRODUCTION 



Split Supersymmetry (SS) was introduced in order to avoid some of the most notorious inconveniences 
of the Minimal Supersymmetric Standard Model (MSSM), namely, the lack of an automatic mechanism to 
avoid large flavour changing neutral currents and CP violation, and fast proton decay Uj]. The strategy is 
to consider all scalars, with the exception of one Standard Model (SM) like Higgs boson, with a very large 
mass of the order of m, using unification of gauge couplings and the lightest supersymmetric particle (LSP) 
as Dark Matter candidate, as the only guiding principles ||2|]. In addition, motivated by the Cosmological 
Constant fine tuning problem, the electroweak scale fine tuning is accepted as a property of nature to be 
explained later by other principles to be discovered. 

In this SS scenario, the light supersymmetric Higgs boson will have SM-like couplings and will be 
difficult to differentiate the two models in the absence of other signals ij]. The Large Hadron Collider will 
shortly start accelerating protons, and the best chance in this case for the larger detectors ATLAS l^j] and 
CMS 1 5] to detect supersymmetry is in the chargino and neutraUno sector |^. While the lightest neutralino 
is the LSP, which is stable and candidate to dark matter in this R-Parity conserving model, the heavier 
neutralinos will decay into it. In SS, X2 ^'^^ have the chance to decay via intermediate scalars, and 
it will do it via intermediate Z bosons, X2 — * Xi^* — ^ Xiff- The other important decay mode is 
generated at one-loop, the radiative decay of the neutralino X2 — * Xi7' where all virtual charged particles 
contribute inside the loop This decay mode is well studied in the MSSM, and despite being generated 
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at one-loop, it can lead to large branching ratios ||8|]. 

In this article we are interested in the one-loop two-body decay mode X2 — ^ Xi7 in Split Supersym- 
metry, and its relative size with respect to the tree-level three-body decay mode X2 — ^ Xiff- We study the 
region of parameter space where the radiative decay is enhanced, showing it to coincide with a relatively 
wide strip around M2 ~ Mi . The signal for the radiative decay, an energetic photon plus missing energy, 
is clean and experimentally attractive, as long as the photon does not become too soft due to lack of phase 
space. We show that in this strip of parameter space, where the photon is still easily detectable 1I4], a mea- 
surement of the two main branching ratios can give information on the supersymmetric scale fh. This is not 
a small feature because it is very difficult to measure the split supersymmetric scale in these models. 

II. SPLIT SUPERSYMMETRY 

Above the scale m the supersymmetric lagrangian is governed by the following R-Parity conserving 
superpotential, 

Wmssm = -X'^HuQU + X'^HdQD + X'HdLE - fiHuHd (1) 

where A", A"', and A*^ are the Yukawa coupling 3x3 matrices, and fi is the Higgs supersymmetric mass 
parameter. ^From this superpotential we highlight the following terms in the MSSM lagrangian, 

Cmssm = -^{Hia^Hu + HyH^f - ^ [hIHu - H^Hdy 

~iml^ + ti^)HiHu - {ml^ + I^'WIh^ + mljH^eHu + h.c.) 
+XlHleu,qj - XfjHjed^qj - X^^Hjeedj (2) 



-^{(ja^W' + g'B) Hu-^ {ga'^W + g'B)Hd + h.c. 



where e = ia2- This lagrangian is valid at the scale fh and above. At fh the Higgs potential is character- 
ized by quadratic terms proportional to squared gauge coupling constants, g{fh) and g'{fh), plus three mass 
terms. The two Higgs eigenstates are found to be rotations of Hu and Hd by an angle (5, the hghtest one 
given hy H = — cos /3ei7^ + sin /3i^„- Also in the MSSM lagrangian we have the Yukawa interactions with 
couplings A"j(m), Xfj{fh), and Xfj{fh). Finally, we see that the higgsino-Higgs-gaugino vertex are propor- 
tional to the gauge couplings, as Higgs-Higgs-gauge boson couplings are, as dictated by supersymmetry. 
The Split supersymmetric lagrangian, valid at the scale fh and below, is given by, 



hfjqjUieH* + h^jdiH + hldjCiH + 



+ ^ [gui'^W'' + g'uB)Hu + ^ ( " dd^y'^W'' + gdB)Hd + h.c. 



(3) 
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where the Higgs field H is the surviving Higgs doublet at low energies. The Higgs potential is defined by a 
mass term m? and a quartic self coupling A. The electroweak symmetry breaking occurs since > 0, and 
the Higgs field acquires a vacuum expectation value v. The matching condition for the Higgs self interaction 
at the split supersymmetiic scale m is 



\{m) 



g^{m)+g'^{m) _2 



cos^ 2(3, 



(4) 



and this coupling should be run down to the weak scale to find the correct electroweak symmetry breaking. 
The Yukawa couplings in the split sumersymmetric model are called hfj , hf- , and hf^ , and at the scale 



the coiTcsponding matching condition are 



h^Am) = A^(m)sin/3, 



Kfi'^) = A-:.'^(m)cos/3. 



m 



(5) 



Finally, we notice from the Split Supersymmetric lagrangian in eq. ^ the Higgs-gaugino-higgsino inter- 
actions, whose couplings have the following matching conditions with the analogous terms in the MSSM 
lagrangian of eq. ©, 



9uim) = g{m) sin P, 
a'ui^) = g' {fn) sin P, 



Qdim) = g{m) cos P, 
a'di^) = 9'{fh) cos/?, 



(6) 



The renormalization group equations for these and other couplings can be found in ref. ||2|]. For illus- 
tration we show in Fig. [T] some of their behaviour. In Fig [Ua) we have the running of the gauge coupling 
constants in Split Supersymmetry, which unify at the GUT scale as in the MSSM. In the second frame 
Fig. [Hb) we plot the Higgs-gaugino-higgsino couplings as a function of tan P for rh = 10^^ GeV. If 
tan P = 1, neither the boundary condition nor the ROE differentiate between g„ and g^ or between g!^ 
and g'^. A sharp splitting appears when tan/3 increases. The down couplings become smaller than 0.1 for 
tan P ^ 10, while the up couplings approach asymptotically a maximum as tan P increases. In Fig. [He) we 
plot the Higgs-gaugino-higgsino couplings as a function of the split supersymmetric scale m. In Fig. (Hd) 
we plot the Higgs-gaugino-higgsino couplings normalized by the gauge couplings evaluated at the weak 
scale as a function of the scale m. We choose the value tan P = 50, and observe deviations up to ±20%. 
Of course, if the split supersymmetric scale is taken equal to the weak scale, there is no deviation. 

Now we introduce the following notation, 

cos 2/3 



tan P = 



tan/3' = 4 
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9d 

9'd 



tan P 



mw 



tan P 



1 + 
1 



647r2 

cos 2/3 
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FIG. 1: (a) Gauge coupling unification is preserved in Split-SUSY. (b) Split-Susy couplings dependence on tan/3, 
(c) Split-Susy couplings evaluated at weak scale as function of m. (d) Split-Susy couplings normalized with gauge 
couplings evaluated at weak scale, for different values of m. 



where it is understood that /3 is defined at the scale m, while f3 and /?' are defined at the weak scale. The 
approximated expressions in eq. © is obtained from the corresponding RGE. These definitions together 
with, 

9^ = alii^w) +9d{^w), g''^ = g'u{mw) +g'/{mw), (9) 

allow us to write the neutralino and chargino mass matrices in such a way it resembles those of the MSSM. 
The mixing angles f3 and f3' are plotted in Fig. |2]as a function of fh for tan /? = 10, 50. Of course, there is 
no difference between the three angles if m = M^^ak- 

With this notation, the neutralino mass matrix in Split Supersymmetric models has the following form. 





Ml 
















M2 


\gvcp 


-\gvsp 


(10) 


-Ig'vc'p 


\gvcp 









\g'v~s'p 


-^gvsp 










which is written in the usual basis = {B,W^,H^,H'^). We have used the notation = cos 
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FIG. 2: Dependence of tan /? and tan (3' on the Split Supersymmetric scale, for different values of tan (3. 

cp = cos /3, s'p = sin/3', sp = sin/3. This mass matrix is diagonalized by the matrix N, such that 
iV*M^,f A^"^ = {M.^^)diag, and the eigenvectors x° = Nijil^^- are the neutrahnos. 

ni. NEUTRALINO DECAYS IN MSUGRA AND SPLIT SUPERSYMMETRY 




FIG. 3: Tree-level three body decay diagrams in Split-SUSY for the second lightest neutralino. 



In Split Supersymmetry the three body decay modes of the second lightest neutralino receive contri- 
butions from intermediate gauge and light Higgs bosons, with negligible contribution from sfermions and 
heavy Higgs boson. These graphs are in Fig.|3l where the major contribution is from the Z-boson exchange, 
since the fermions in the final states have a very small mass. We calculate these decay rates integrating over 
the phase space with numerical techniques. We compare our calculations for the case of a small SS scale 
m ~ 1 TeV with results from the ISASUGRA code ||9|] with mo ^ 1 TeV. These are in agreement within 
small differences, the main of which is the distinctive running of Higgs-higgsino-gaugino couplings present 
is SS. 

Our main interest in this article is the radiative decay of the second lightest neutralino into the LSP and 
a photon. This one-loop generated decay can shed light into the properties of heavy particles present in the 
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loop, charginos in the case of Split Supersymmetry. In addition, it is an experimentally interesting decay 
mode since includes only a hard photon plus missing energy. Contributing loops in SS are displayed in 
Fig. m The loops include both charginos and W gauge bosons (charged Goldstone bosons are implicit). 
AU other charged scalars which could contribute have a mass of the order of m and they are neglected. Of 
course, the effect of the heavy particles is felt via the RGB of the effective couplings below m. We calculate 
the integral over internal momenta analytically using dilogarithms i^. 
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FIG. 5: Branching ratio and decay width for the different neutralino decay modes, in mSUGRA as a function of mo 
and around the SPSla benchmark point. 



In Fig. [5] we show the behaviour of the three main decay modes of the second Ughtest neutralino in 
mSUGRA: the tree-level two-body decay X2 ~^ f-^i^, the tree-level three-body decay X2 ~^ Xiff' the 
one-loop two-body decay X2 ~^ Xi7- have Branching Ratios to the left and decay rates to the right, as 
a function of the universal scalar mass mQ. The tree-level modes are calculated with the code ISASUGRA, 



while the one-loop mode with our code. The other parameters are taken as in benchmark SPSla 1110(1 . which 
is given in Table Jl 
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TABLE I: Input parameters for SPS 1 a mSUGRA benchmark point. 



Parameter Value Units 



7710 


100 


GeV 


Ml/2 


250 


GeV 


Ao 


100 


GeV 


tan (3 


10 




sign(At) 


+ 1 





TABLE II: Chargino and Neutralino masses for SPS la. 



Particle Mass Units 


X? 


99 


GeV 


xl 


175 


GeV 




352 


GeV 


xl 


372 


GeV 


xt 


175 


GeV 


xt 


372 


GeV 



We can see that the one-loop decay rate r(x2 Xi7) has a maximum near tuq = 150 GeV, decreasing 
for larger universal scalar mass. At values above ttiq = 1 TeV the mSUGRA result differs between 3- 
20% from our own SS result, calculated with values m ^ 1 TeV. Of course, it is expected that mSUGRA 
results approach the SS result for large mo, since at large values of the universal scalar mass the triangular 
contributions from heavy scalar particles diminish. In SS the effect of these heavy particles appear through 
the RGE of the different couplings, but small differences remain between the two approaches because the 
running couplings include leading logarithm effects from all loops. We remind the reader that these RGE 
effects in SS do not spoil the unification of gauge coupling constants, as stressed in ref. and illustrated 
in Fig. [T] The branching ratio B{x2 ~^ Xi7) remains between 10"^ and 10"^ for this mSUGRA scenario. 
In Table |II] we show the neutralino and chargino masses for ttiq = 100 GeV, with masses only slightly 
increasing (one or two GeV) for larger scalar mass, calculated using SUSPECT inl]. 

We compare in Fig. [5] the one-loop generated decay mode X2 ~^ Xi7' with the tree-level decay 
modes. We define the tree-level two-body decay X2 ~^ (-^t^ as the sum of the three leptonic decays 
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TABLE III: Slepton and squark masses for SPSla. 



Particle Mass Units 
eR,fiR 145 GeV 
eL,p-L 204 GeV 





136 


GeV 




208 


GeV 


ii 


375 


GeV 


bi 


491 


GeV 



TABLE IV: Input parameters for Split Supersymmetry benchmark point. 



Parameter Value Units 



Ml 


102 


GeV 


Ma 


192 


GeV 


Ms 


587 


GeV 


M 


357 


GeV 


tan (3 


10 





X2 ee, /i/x, fr which occurs for values mo ^ 160 GeV, where the sleptons have a mass smaller than 
m^o. In this region, this decay mode dominates with a branching ratio near unity. We also have the tree- 
level three-body decay X2 ~^ X?//' where we sum over all possible fermions. Above mo ~ 160 GeV 
the off-shell intermediate particles which contribute are the Z gauge boson and the squarks or sleptons, 
depending whether the final state fermion is a quark or a lepton. In this region the B{x2 Xiff) is near 
unity. Below mo ~ 160 GeV the contribution from the intermediate light on-shell sfermion is removed, and 
B{x2 ~^ Xif f) drops to a value between 10"'^ and 10~^. In Table Hill we show the slepton and the lightest 
squark masses for the SPSla scenario using SUSPECT. For larger mo these masses grow up sharply, with 
right selectron and smuon becoming on-shell if mo ^ 150 GeV and similarly for the lightest stau if mo ^ 155 
GeV. Both thresholds are fused into one in Fig. [5]because of the low resolution used in the graph. 

In Fig. |6]we show the tree-level three-body decay X2 ~^ Xiff ^nd the one-loop two-body decay X2 ~^ 
Xi7 in Split Supersymmetry, as a function of the wino mass M2, both calculated with our own code. We 
choose as SS benchmark point the one given in Table HVl whose soft gaugino and higgsino mass values 
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Ma (GeV) Mg (GeV) 



FIG. 6: Branching ratio (a) and decay width (b) for X2 decay modes in SS as a function of M2. In (c) we see the 
mass eigenvalues crossing between the second and third neutralinos, which cause the discontinuous behaviour of the 
branching ratio and decay width (d). 

coincide with the low energy soft masses from SPSla. When varying the wino mass M2 in Fig.[6j we vary 
also the bino mass Mi keeping constant the M2/M1 ratio as in our SS benchmark scenario. In frame (a) we 
have the branching ratios, where tree-level thi^ee-body decay dominates over the one-loop decay X2 ~^ Xil 
with a BR near unity. Since we work in SS, the X2 ~^ Xif f mode is mediated only by an intermediate 
Z gauge boson. Similarly, the X2 ~^ Xi7 decay is generated only by quantum corrections with W gauge 
bosons and charginos inside the loop. 

In frame (b) we plot the decay rates for these two modes. The discontinuity on both decay rates occur 
near M2 ~ 560 GeV and corresponds to an eigenvalue crossing. In frame (c) we see this eigenvalue 
crossing, with a X3 higgsino type and X2 gaugino type for M2^560, while the opposite occurs for M2^560. 
The effect of the crossing can be seen very clearly in frame (d) where we have the three-body decays for 
both neutralinos X2 ^^'^ Xs- 
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IV. HIGH RADIATIVE NEUTRALINO DECAY BRANCHING RATIO 



In this chapter we analyze with more detail the radiative decay for the second lightest neutralino, and 
look for conditions for an enhanced B{x2 ~^ Xil) in SpUt Supersymmetry. In Fig. |7^ we show with a 
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FIG. 7: (a) Radiative decay width, (b) 3-Body decay width, and (c) Branching ratio for the radiative mode in M1-M2 
plane in Split Supersymmetry. In (d) we have different curves as a function of M2 showing the special behaviour at 

Ml = Af2. 



color code the logarithmic values for the decay rate T{x2 ~^ Xil) ii^ the gaugino mass plane M1-M2, with 
r measured in GeV. We show gaugino masses smaller than 500 GeV, and vary randomly the values for 
and tan (3. The largest values for the decay rate occur for Mi > M2, reaching typically up to 10^^ GeV. In 
the opposite case, when Mi < M2, the decay ratio varies typically between 10"^ and 10"*^ GeV. There is a 
narrow fissure around Mi ~ M2 where the decay rate drops to values between 10^^*^ and 10^^^ GeV. The 
fissure is not exactly at the bisector but somewhat deviated to the Mi > M2 side of the quadrant. In this 
fissure Am^ = m^o — m^o is minimum, and the drop of the decay rate is a kinematical suppression. 

In Fig. |7j5 we have a similar plot for the decay rate r(x2 Xif f)^ which is much larger than the 
previous case. This decay rate is more symmetrical with respect to the bisector, with decay rates reaching 
maximum values between 10"^ and I GeV, and a deep fissure at Mi ~ M2 going all the way down to 10~^^ 
GeV or smaller. The fissure is situated over the bisector, and it is due to a zero in the neutralino-neutralino-Z 
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coupling, i.e. a dynamical suppression. 

The branching ratio for the radiative decay B{x2 ~^ Xil) is shown in Fig.|7]:. We see that the second 
lightest neutrahno one-loop generated decay X2 ~^ Xil can dominate in a wide zone near the bisector 
Ml ~ M2. The reason for this possibility is that both decay rates decrease shaiply in parallel fractures but 
slightly displaced. One of the fractures due to a zero in the x^X^^ coupling (dynamical), and the other due 
to an eigenvalue degeneracy (kinematical). This is confirmed in Fig. |7]l, where we show both decay rates, 
the radiative decay branching ratio, and the x^X^^ coupling as a function of M2, with constant values for 
Ml = 275 GeV, ^ = 400 GeV, tan (3 = 10, and m = 10"^ GeV. We see that the zero for T{x^ XifJ) 
coincides with the zero for x^X^^ coupling at Mi = M2, and that the minimum for r(x2 XiT) coincides 
with the point where m^o — m^o is minimum, at A/2 ^ Mi. 

In order to better understand the above result, it is instructive to see the neutralino mass matrix in the 
basis [—ij, —iZ, Hi,H2], where 



(11) 



7 




cw sw 




B 




Hi 




0/3 


-S/3 




Hd 


Z 




—sw Cw 




W 




H2_ 






C/3 _ 




Hu 



In this basis the mass matrix in eq. (flOl ) looks as follows. 



(12) 



with M^o = (M^o)^- The different submatrices are equal to. 



Mic^ + Mas^ (M2-M1 

)swcw 

{M2-Mi)swcw Mislr + M2c'^ 



flS2/3 -/UC2/3 
-HC2f3 -/iS2/3 



(13) 



for the blocks in the diagonal, and 



(14) 



Cjsigdsw - g'dCw) + s/sigusw - g'ucw) spigdsw - g'dCw) - ci3{gusw - g'uCw) 
cwicpgd + spgu) + sw{cpg'^ + s0^) cwispgd - cpgu) + swispg'd - cpg'u) 

for the off diagonal block. This neutralino mass matrix reduces to its analogous expression in the MSSM if 
we neglect the running from m and the weak scale: 





g/cw 



as 



(15) 
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In the MSSM case, the direct mixing between photino and higgsinos vanishes, but a direct coupling be- 
tween zino and one of the iiiggsinos remains. This implies that in general the lightest neutralino has a non 
vanishing component of higgsino, which in turn translates into a non vanishing x^X^^ coupling. In this 
way, the photino will decouple from higgsinos in the region Mi M2, as seen from eq. ([T3l) . and the decay 
rate T{x2 ~^ Xiff) vanishes also, making the decay X2 ~^ Xi7 the dominant one. In SS the mechanism is 
similar, but modified by RGE effects. 

As we mentioned, in the dynamical suppression region where Mi ~ M2 the decay mode r(x2 —>■ Xiff) 
is suppressed because the Z coupling to the photino is absent, thus, when the LSP is nearly photino, the 
ZX1X2 coupling is neai^ly zero. In the MSSM this region does not exactly coincides with the kinematical 
suppression region where m^o m^o due to the remanent higgsino-gaugino mixing seen in eq. (fT5l) . In 
this case, phase space is small, and X2 he forced to decay into Mght mesons. The situation is similar 

>■ 10-6 

10-9 
10-12 

10-15 



20 25 30 35 40 45 50 

Am^ (GeV) 

FIG. 8: Decay rate for — * Xi'''"*' in comparison to the decay rate for X2 — * Xiff^ a function of the mass 
difference Am^ = — rriy,-^ . 

in Split Supersymmetry where the difference lies in the fact that in SS RGE effects separate further the 
regions where ZX1X2 coupling and Am^ vanish, as indicated by the higgsino-gaugino mixing in eq. (fT4b . 
In Fig. [8] we have plotted the decay rate T{x2 — ^ Xiff) ^rid the (included in the former) decay rate 
r(x2 — ^ Xi^*^) a function of Am^. The independent variable we are varying is M2, exactly as in 
Fig. |7}i. Both decay rates vanish at the point where the ZX1X2 coupling is null, but at that point A.m^ is 
not zero. Indeed, in Fig. ^ we have Am^ = 26.8 GeV for m = 10 TeV, while RGE effects changes it to 
Am^^ = 29.2 GeV for fh = Mqut- Therefore, the photon in the decay X2 ~^ Xil will have enough energy 
to be easily detected. Note that the two branches in each decay are defined by the sign of M2 — Mi . 
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As we discussed, in Split Supersymmetry the mechanism is analogous to the MSSM, but the details are 
modified by the Renormalization Group Equations effects. Indeed, the remaining higgsino component of 
the lightest neutralino in the case Mi = M2 is controlled by the SS scale fh via the RGE effects on the 
different couplings. This can be seen in Fig. |9] where we have the X2 branching ratios dependence on tan /3 




FIG. 9: X2 decay modes as a function of the Split Susy scale m and tanj] for the scenario where A/i = M2- 



and the SS scale rh, with /i = 400 GeV and Mi = M2 = 275 GeV. The dependence on tan /3 is relatively 
mild, as it is the dependence on fh. But in comparison to other observables, the dependence on the SS scale 
is very important, because in this scenario a measurement of the X2 branching ratios could yield valuable 
information on the SS scale, otherwise difficult to extract from experiments. In the right frame we see that 
B{y^ x\l) dominates over a large region, with a small decrease at small tan (3 and large m. In the left 
frame we have B{x2 ~^ Xiff)^ with values that go from 10~^ up to 0.2. Clearly, a measurement of the 
branching ratio can give valuable information on the split supersymmetric scale. 
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V. CONCLUSION 

We have calculated the decay rates and branching ratios of the second lightest neutralino in Split Su- 
persymmetry, and concentrate in the radiative decay X2 ~^ Xi7- We compared our results with mSUGRA, 
finding agreement when the scalar mass parameter is very large, mo ~ 1 TeV, where B{x2 ~^ Xil) 10~^. 
Small differences remain due to RGE effects, which increase with larger split supersymmetric scale m. For 
larger values of mo comparison is not possible since large squark masses in quantum corrections destabilize 
the EWSB in direct mSUGRA calculations. 

In general models, the possibility that m^o is not much different than the mass of the LSP is experimen- 
tally challenging. This is because the decay products that can be detected, photons in the case of X2 ~^ Xi7 
and fermions in the case of X2 ~^ soft and specialized analysis have to be done with the data. 

Nevertheless, in our model the radiative decay dominates in a region where Am^ is large enough to produce 
an energetic photon. We focus on the region Mi ~ M2, where this possibility is realized and show that the 
decay X2 ~^ Xi7 is dominant in a wide band around the bisector Mi = M2. Furthermore, in this region a 
measurement of the branching ratios B{x2 ~^ Xil) ^rid B{x2 ~^ Xiff) can give information on the value 
of the split supersymmetric scale m. 
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